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ABSTRACT 



NGC 6251 is a luminous radio galaxy as 104 Mpc away that was detected significantly with the Fermi Gamma-ray Space Telescope, 
and before that with EGRET (onboard the Compton Gamma-ray Observatory). Different observational constraints favor a nuclear 
origin for the y-ray emission. Here we present a study of the spectral energy distribution (SED) of the core of NGC 6251, and give 
results of modeling in the one-zone synchrotron/SSC framework. The SSC model provides a good description of the radio to y-ray 
emission but, as for other misaligned sources, predicts a lower Lorentz factor (T ~ 2.4) than typically found when modeling blazars. 
If the blazar unification scenario is correct, this seems to point to the presence of at least two emitting regions in these objects, one 
with a higher and one with a lower Lorentz factor. The solution of a structured jet, with a fast moving spine surrounded by a slow 
layer, is explored and the consequences of the two models for the jet energetics and evolution are discussed. 

Key words, y-rays: galaxies — galaxies: active — galaxies: jets — galaxies: individual: NGC 6251 



1. Introduction 

The Fermi Gamma-ray Space Telescope has confirmed mis- 
aligned AGNs (MAGNs), including radio galaxies (RGs) and 
steep-spectrum radio q uasars (SSRQs), as a new and important 
class of y-ray emitters (lAbdo et al.ll2010cfej). In the fir st year of 
activity, the Large Area Telescope (LAT, Atwood et al. 2009) on- 
board Fermi detected 1 1 MAGNs belonging to the 3CRR, 3CR, 
and MS4 catalogs at 178 and 408 MHz. At these low frequen- 
cies, the Cambridge and Molonglo surveys favor the detection of 
radio galaxies over blazars. Seven of the 1 1 Fermi MAGNs are 
nearb y (z ^ 0.06) Fanaroff-Riley type I (FRI; Fanar off & Rilevl 
' 1 19741) radio galaxies, while the rest are FRII radio galaxies. 
The FRI radio galaxies detected with Fermi are significantly less 
y-ray luminous than BL Lac objects. The BL Lac objects have 
isotropic 100 MeV - 10 GeV luminosities L y * 10 44 - 10 46 
erg s , significantly higher than the FRI radio galaxies, which 



have L y 



10 4 



10 44 ergs-'. 



They also have spectral slopes 



that are consistent with low- or intermediate-synchrotron-peaked 
blazars (Abdo et al. 2010c a). This is consistent with AGN uni- 
fied schemes, according to which an increase in the inclination 
angle of the jet axis with respect to the observer's line-of-sight 
implies a de-amplification of the observ ed flux and a general 
shift in th e emission to lower energies (lUrrv & Shaferl fl984: 
Urry & Padovani 1995, and references therein). The classifica- 
tion itself of misaligned sources is relative to our expectations 
of detecting in radio galaxies strongly beamed jet radiation from 



highly relativistic (F > 10) bulk flows, as in the case of blazars. 
Thus, for the sake of clarity, by misaligned here we do not mean 
that the observer line-of-sight is outside the radiation cone of the 
source but rather that it falls outside the (narrow) beaming cone 
we would expect in the case of a highly relativistic core. 
A full comparison with the expectations of the unification sce- 
narios for AGNs has not yet been considered, but the addition of 
a substantial number of misaligned AGNs at 100 MeV - 10 GeV 
energies now enables this to be made. 

Here we present our study of the SED of the core of NGC 
6251, identifie d with 1FGLJ1635 4+8228 in the First Fermi - 
LAT Catalog (lAbdoetal.1 l2010bllcl First Fermi -LAT Source 
Catalog, 1FGL, and First Fermi -LAT AGN Catalog, 1LAC, re- 
spectively). NGC 6251 is the fifth y-ray brightest radio galaxy 
in the MAGN sample with an integrated flux F_ 9 = (36 + 8) 
in units of 10~ 9 ph(> 100 MeV) cirT 2 s" 1 . It is associated wit h 
the EGRET source 3EG J 162 1+8203 (lMukherieeetal.1120021) . 
being after Centaurus A one of the most likely associations 
of an EGRET source with a radio galaxy (also of interest are 
NGC 1275, 3 C 111, and M87). B eing bright and relatively close 
(z = 0.0247 l lWegner et al]l2003l NGC 6251 provides a nearby 
cosmic laboratory to explore the jet structure of radio sources 
hosting a supermassive black hole. 

After describing the source in more detail (Sect. 2), we describe 
the multiwavelength datasets used to assemble the spectrum of 
NGC 6251 (Sect. 3). In Sect. 4, we illustrate some uncertainties 
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regarding jet orientation and location, before presenting the re- 
sults of one-zone synchrotron self-Compton (SSC) modeling in 
Sect. 5. There we derive magnetic fields, outflow Lorentz fac- 
tors, and absolute power estimates for comparison with blazars. 
Parameter comparisons between SSC models of BL Lacs and 
FRI galaxies reveal the incongruity between these two classes 
of sources, which we discuss in light of the unification scenario 
and explore the feasibility of an alternative scenario. In Sect. 6, 
we consider the connection between the jet power, the disk lu- 
minosity, and the accretion. The results are summarized in Sect. 
7. 



2. NGC 6251 



NGC 6251 is classified as a FRI radio galaxy dLaing et al.ll 1983b 
based on of its morphology and monochromatic radi o luminos- 
ity of * 0.94 x 10 31 erg s" 1 sr" 1 at 178 MHz (Wagg ett et alJ 
1977)Q. Nevertheless, it exhibits an unusual radio morphology, 
with an sa 1 ° extensi on of bright radio emission to the northwest 
dWaggett et al.ll977b . and characteristics typical of both FRI and 
FRII radio galaxi es (see the radio maps at 327 and 608 MHz in 
Mac k etal.ll 1997b . At a distance of a 104 Mpc, this corresponds 
to a linear extension of « 1 .8 Mpc (1" = 0.5 kpc). 
High resolution 1.48 and 4.9 GHz VLA studies reveal a radio 
core and a complex jet, that is bright and structured within 4.4' 
of the core (~ 113 kpc), and then faint and curved at larger 
size scales (Perlev et al. [|1984b . On the basis of radio images, 
Sudou & Taniguchi (2000) showed that the angle between the 
jet axis and the line of sight increases by ^ 10°, going from 
33° at 50" from the core to 45° at 200". The VLBI maps show 
an asymmetric core-jet radio structure aligned with the VLA jet 
(Jones et al. 1986). 

Both non-thermal and thermal SED components contribute to 
the core emission of NGC 6251. At wavelengths between * 15 
and 30 /mi, thermal dust emission appears to dominate over 
an estimated ~30% nonthermal contribution to the total mid-IR 
Spitzer flux (Leipski et al. 2009). Synchrotron radiation proba- 
bly also accounts for the bulk of the optical to UV emissio n 
within a 0.2" of NGC 625 l's core dChiaberge et all [2003b . 
The nuclear region in the optical band is complex. A warped 
dusty disk « 1.43" in extension unevenly reflects UV radia- 
tion f rom the nucleus, as see n in Hubble Space Telescope im- 
ages (Ferrarese & Ford 1999). Ionized gas in the 0.3" region 
(a 150 pc) surrounding the nucleus implies that the nucleus of 
NGC 6251 harbors a bl ack hole with mass * (4 - 8) x 10 8 M 
dFerrarese&Fordlll999h . 

NGC 625 1 has been observed with space-based observatories at 
X-ray energies. The unresolved core is the main X-ray emitter, 
although high-resolution imaging with Chandra resolves distinct 
X-ray emission in three different jet regions (Evans et al. 2005). 
The detection of the Fe Kg line is still debated (Turner et al 



1997[ iGuainazzi et all 120031: iGliozzi et alJ l2004t lEvans et al 
2005b . An extended, ~ 100 kpc X -ray (k B r a 1.7 keV) thermal 



halo was detected with ROSAT (Birkinshaw & Worrall 1993; 
Mac k et al.l Il997h. Chandra (lEvans et al. 2005b, and XMM- 
Newton (Sambruna et alJl2004b . A drop in the surface bright- 
ness of the X-ray halo, in positional agreement with the northern 
radio lobe, suggests that the lobe has evacuated a cavity in the 



1 Throughout the paper, we adopt the cosmological parameter values 
H = 71 km s- 1 Mpc" 1 , Q ffl = 0.27, and n A = 0.73. 

2 A ccording to the classification proposed by Fanaroff & Riley 
J19741) . the flux density at the FRI/FRII boundary, corrected for the 
adopted cosmological values, is 1 x 10 32 erg s sr -1 . 



surrounding X-ray gas dKerp & Macldl2003t lEvans et al.ll2005b . 
NGC 6251 has also been detected at hard (> 10 keV) X- 
ray energies by Beppo-Sax (Guainazzi et al. 2003; Grandi et al. 
2006) and INTEGRAL dFoschini et al.ll2005l) and, as previously 
noted, was propos ed as a counterpart to the EGRET source 
3EG J1621+8203 dMukherjee et all [2002b . Obser vation from 
Fermi confirmed that NGC 6251 is a GeV source (Abd oetall 
20101©. The origin of the nuclear SEP of NGC 6251 
has been discu s sed in several papers (O iiaberge et al. 2003; 
Guai nazzi etaP 12003b iGhisellini etail 120051: iFoschini et"atl 
2005), which have concluded overall that the observed nuclear 
emission is likely to be dominated by emission from a relativis- 
tic jet. The observed SED indeed shows the typical double-hump 
shape characteristic of blazar sources. 



3. SED Data 

3.1. Radio to optical/UV data 

The radio to UV data collected from literature and shown in 
Table |4] are used to assemble the nuclear broadband SED of 
NGC 6251. We note that radio fluxes are taken at different 
angular resolution and sample nuclear regions of different 
dimensions. Moreover, the data are not simultaneous and 
variability could be an issue. lEvans et all (|2005), however, 
inspected VLBI core radio fluxes over 17 years and limited the 
maximum radio fl ux variability to a factor of < 2 (see Table 9 in 
lEvans et al.l2005l) . Data in the microwave region are provided by 
the five-year Wilkins on Mic rowave Anisotropy Probe (WMAP) 
catalog (Wrig ht et al.ll2009h and by the published Planck Earl y 
Release Compact Source list (Planck Collaboration et al. 201 1). 
The IR flux (15-30 pm) is decomposed into ~30% synchrotron 
jet emission extrapolated from the radio data and the remainder 
thermal radiation (Leipski et al. 2009). 

The non-thermal origin of the bulk of the nuclear optical 
emission is supported by the high degree of polarization of the 
UV emission (close to 50%) and the position of NGC 6251 
in diagnostic diagrams for radio galaxies, in common with 
the other FRI sources (Chiaberge et al. 2003, and references 
therein). The drop in the optical-UV band might be either real or 
cause by the dust extinction of the intrinsic optical and UV flux. 
It is important to establish wether the IR to X-ray flux can be 
ascribed to a single non-thermal proces s. The interpretation of 
the available data is partly controversial. IChiaberge et all {2003) 
argued against significant dust reddening because the progres- 
sive steepening of the optical-UV slope, related to the increase 
in the dust extinction with frequency, is not observed. By visual 
inspection, the optical nucle us appears to be very bright and 
unobscured by the dust lane. lFerrarese & Fordl (11999) estimated 
a mean intrinsic visual extinction (Ay) of 0.61 ± 0.12 mag. 
The total extinction, which also accounts for the foreground 
reddening in the direction of NGC 6251, is A v - 0.88 ±0.13 
mag. This value indeed relies on some assumptions, such 
as that the reddening law follows Card elli et al.l ([1989) and 
R v = A V /E(B - V) is equal to 3.1. An upper limit of 1.0-1.5 to 
the value of Ay was derived in lChiaberge et al.l d2003h . 
In the SED (Figure @}, we show both the optical fluxes for no 
absorption (black solid triangles) and after de-reddening (empty 
circles) for Ay = 0.88 mag (using the extinction curves of 
Cardelli et al. 1989). The SED has a clear double-peaked shape. 
However, even in the case of absorption (Ay < 1.01), a unique 
emission process for the IR to X-ray data seems unlikely as the 
IR and optical fluxes do not lie on the extrapolation o f the X-ray 
slope to the lower frequencies (see also Figure 13 in lEvans et al.l 
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I2005h . 



3.2. High energy data 

All data from XMM-Newton, Swift, and Chandra observations 
available in the public archives were re-analyzed and included in 
the SED. For EGR ET and Fermi-L AT data, we use the sp ectra 
provided by Mukh eriee etall (120021) and lAbdoetail (1201 Pel) . 



3.2.1. X-ray spectral analysis 

The XMM-Newton observation of NGC 6251, performed in 
March 2002, was analyzed using the SAS v.9.0 software and 
available calibration files. We excluded time intervals affected 
by high background. After this data cleaning, we obtained a net 
exposure and count rate of 8.7 ks and 1.867+0.015 count s 
for the pn, 13.9 ks and 0.553+0.006 count s _1 for MOS1, and 
13.5 ks and 0.556+0.006 count s" 1 for MOS2, respectively. The 
source and background spectra were extracted from circular re- 
gions of 27" radius. The response matrices were created using 
the SAS commands RMFGEN and ARFGEN. The nuclear data are 
not piled up. Data were grouped into 25 counts per bins in order 
to be able to apply the x 2 statistic. The best-fit model for the 
pn (0.3-10 keV) data consists of an absorbed power-law plus 
an APEC component which models a collisional gas emission 
(^ 2 =410 for 372 degrees of freedom). The parameter va l ues, re - 
ported in Table [1] agree with the results of [Evans et al.l ([2005). 
The same results (not shown in Table [TJ were obtained using 
MOS data. There is no significant evidence of the Fe Ka emis- 
sion line in the XMM data. 

The Swift X-Ray Telescope observed NGC 625 1 three times be- 
tween April 2007 and May 2009. The X-ray data were reduced 
using the on-line XRT data analysis provided by the ASDQj. 
Source spectra for each observation were extracted from a circu- 
lar region of 20" radius, while the background was taken from an 
annulus with an inner radius of 40" and outer radius of 80". The 
data were re-binned to 20 counts per bin in order to calculate^ 2 . 
All spectral fits were performed in the 0.5 - 10 keV band. Data 
can be well described by an absorbed power-law with column 
density slightly in excess of the Galactic value {Nh Gi ,, - 5.4xl0 20 
cm -2 ). Unlike XMM Newton, the XRT spectra, characterized by 
a lower signal-to-noise ratio, do not require the addition of a soft 
thermal emission component (see TableQ]). We note that, during 
the 2002 observation, the source appeared in a higher state (see 
Sect. 3.2.2) 

We analyzed only the most recent and longest Chandra obser- 
vation of the source in November 2003 that was performed in 
ACIS-S (S3 and S4 chips) configuration. The Chandra point- 
ing of September 2000 was not taken into account because in 
this observation the core fell in a chip gap. Data were repro- 
cessed using the Chandra Interactive Analysis of Observation 
(CIAO v.4. 1 ) and the Chandra Calibration Database C ALDB 4.1. 
After correction for high particle background, the total exposure 
time is reduced to 45 ks. The nuclear spectrum was extracted 
from a circular region of 5" radius centered on the source. The 
background was chosen in an adjacent circular region with 12" 
radius. Because of the brightness of NGC 6251 at X-ray ener- 
gies, the data are affected by pile-up, which was estimated to 
be w 13% using the PIMMS software! For this reason, we de- 
cided not to consider the Chandra nuclear fluxes in the SED. As 




3 http://swift.asdc.asi.it/ 

4 http://cxc.harvard.edu/toolkit/pimms.jsp 



Energy (keV) 

Fig. 1. Unfolded spectral model of the combined XMM-Newton 
and Swift datasets. Black points: Swift observation of April 
2007; red: Swift observation of May 2009; green: Swift obser- 
vation of June 2009; blue: XMM-Newton observation of March 
2002. 



an extended extranuclear emission associated with the kpc jet is 
clearly visible in the Chandra image, a spectrum was also ex- 
tracted for this component. The accumulation region is a box of 
20" and 5.3" each side, far (« 6") from the nucleus. A box of 
similar size was chosen to estimate of the background. A good 
fit is obtained by an absorbed power law with F=2.3^g 2 4 and 
Nw Co , =5.7x 10 20 cm 2 . The photon index val ue is in good agree- 
ment with the one (r=2.29+° J 4 ,) reported in lEvans et all 12005) 
for the region containing the inner jet (extending from about 1 .5" 
to 27" from the nucleus). The 0.5 - 2 keV unabsorbed flux is 
3.1xl0~ 14 erg enr 2 s _1 , while the 2-10 keV flux is 2.3xl0~ 14 
erg crrT 2 s _1 . 



3.2.2. High energy variability 

We first explored fast (hour-scale) temporal variability of the nu- 
clear X-ray flux during the XMM and Swift observations pro- 
ducing a lightcurve in the 0.5 - 10 keV band for each observa- 
tion. None of the four lightcurves reveal a statistically signifi- 
cant change of flux. To investigate the presence of X-ray flux 
variability on timescales of months/years, we successively per- 
formed a simultaneous fit of the XMM and Swift spectra. The 
four datasets were simultaneously fit with a composite model, 
an absorbed power law, and a thermal (APEC) component, which 
correspond collectively to the best-fit model for the spectrum 
with the highest signal-to-noise ratio (i.e. the XMM-Netwon 
dataset). An energy-independent multiplicative factor (n) allows 
us to compare the fits of the four datasets. The results of the si- 
multaneous fit are reported in the last column of Table [1] The 
values of the model parameters are driven by the XMM dataset. 
As also shown in Fig. [TJ there is an evident flux variation over 
the entire 0.5-10 keV band between the XMM and Swift ob- 
servations, expressed by the different n values. The X-ray flux 
decreases by about 40% between the XMM and the first Swift 
observations and then an additional ~ 15% by the two last Swift 
observations (the values of n parameter for the three observations 
being n XM M = 1, n Sw ifti = 0.60 ± 0.04 and n Sw i/a = 0.45 ± 0.05 
respectively). The close to two years between the first and the 
second Swift observations, performed on April 6, 2007 and May 
5, 2009, respectively, provide us with the shortest interval during 
which X-ray flux variability has been detected. Thus, we con- 
sider this as an upper limit to the X-ray flux variability timescale. 
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Table 1. X-ray data analysis results. 



Satellite 
Obs date 



XMM-Newton/PN Swift/XRTi Swift/XRT 2 Swift/XRT 3 XMM-Swift 
2002-03-26 2007-04-06 2009-05-05 2009-06-05 Combined fit <n) 



k B T (keV) 
normf/' (xl0~ 5 ) 
N„(xl0 21 cnT 2 ) 

r 



normPCxlO" 3 ) 



a {d) 

XMM 



n 



So/if t\ 



n 

" Swift! 



n (rf) 

^(d.O.f) 



0.6±0.2 


- 


- 


- 


0.8^ 


3.0±2.0 


- 


- 


- 


4±2 


0.54+0.01 


98 +0 ° 8 

"• ?o -0.07 


i 4 +o.i 


6 +al 

"■"-0.6 


1.1+0.1 


1.89±0.04 


2.00!°^ 


2 20 +0 - 35 


2 03 +0 35 


1.89+0.04 


1.18+0.05 


0.«5!2S 


071+ 0.27 
u./i I9 


60 +021 
"•""-o.is 


1.18+0.05 
1 


- 


- 


- 


- 


0.60+0.04 


- 


- 


- 


- 


0.45+0.05 


- 


- 


- 


- 


0.44+0.04 


410(372) 


21(29) 


9(11) 


15(14) 


455(430) 



F1ux w 
Flux w 



2.7xl0~ 12 1.9xl0~ 12 1.6xl0~ 12 1.33xl0" 12 

3.6xl0~ 12 2.11xl0 12 1.42X10- 12 1.45xlQ- 12 



Notes. <n) Simultaneous fit of the four observations. <b> Normalization for APEC model is norm = -rrrrr, — ttt f n e n H dV with the abundance table set 

to lAnd ers & Grevessejl[T983). <c) In units of photons cnT 2 s _l . (rf) Energy-independent multiplicative factor for the simultaneous fit. w Unabsorbed 
fluxes for the power-law component in units of erg cirT 2 s~' . 



No time variability in the y-ray band was detected in the 
first 15 months of LAT observations of the Fermi source 
1FG LJ1635 .4+8228, which is associated with NGC 6251 
JAbdo et alJ2010d) . In Fig.EJ we show the 100 MeV - 100 GeV 
lightcurve of the 15 -month LAT dataset generated by dividing 
the total observation period in 5 time intervals of 3-month dura- 
tion. 

We considered longer timescales when comparing the LAT and 
EGRET results. The LAT flux F-g - 36 ± 8 is lower in compar- 
ison with the y-ray flux measured by EGRET, F-9 - 74 + 23. 
However, the large uncertainties in the fluxes make it difficult to 
firmly establish whether this is related to an actual variation in 
the y-ray flux of the source on timescales of years rather than to 
contamination of the EGRET flux from other sources. At a dis- 
tance in the range of ~ 2°-5° from 1FGLJ1635.4+8228, there 
are at least three other y-ray sources in the 1FGL with a flux 
similar to or brighter than 1FGLJ1635. 4+8228, none of which 
has an EGRET counterpart. 



4. Observational constraints 

We must now consider our observational constraints when de- 
veloping a multiwavelength database for modeling, including the 
relative imaging capabilities of different detectors, epochs of dif- 
ferent observations, separate spectral components, and questions 
about whether there is either intrinsic or Galactic extinction. 
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Fig. 2. Fermi -LAT lightcurve of the source 1FGL1635.4+8228, 
associated with the radio galaxy NGC 6251, between 100 MeV 
and 100 GeV. The time is measured from 2008 August 4 and 
covers 15 months of Fermi operation. Each bin corresponds to 3 
months of observations. The horizontal line is the constant flux. 
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4.1. Location of the y-ray emission region 

The construction and interpretation of NGC 625 l's SED de- 
pends on the site of production of the y-ray photons. Because 
of the limited imaging capability of the Fermi-hKY, with a PSF 
of w 0.6° at 1 GeV (Atwood et al. 2009), Fermi can distin- 
guish the central regions from sites coincident with the ex- 
tended radio structures in NGC 6251, but is unable to resolve 
details on length-scales % 3' (1' corresponds to =s 30 kpc). 
Radio galaxies are expected to be sources of y-ray emission pro- 
duced hi_bofli_flie_com£ac_^_core_Mdjivner jet (Marscher et al. 
2008t iFermi-Lat Collaboration et al.ll20loh . as well as emitting 
extended y-ray emission from, at least, photons of the cosmic 
microwave and extragalactic background light (CMB and EBL 
respectively) that are Compton-scattered by the radio-emitting 
electrons in the lobe s (Cheung 2007; Georganopoulo s et al.l 
I2008t lHardcastle et alJl2009h . 

If the source varies, simultaneous multi-frequency VLBI , 
X-ray/Swift, and Fer mi campaigns (S okolovskv et alj 1201 Ot 
iGiovannini et aLlfeOlOh can help to constrain the origin of the 
y-ray emission as well as the physical parameters of the re- 
gion. For blazars, results based on high frequency radio observa- 
tions, rapid y-ray variability, and correlated variability in differ- 
ent wavebands make locations in the vicinity of the bright radio 
core the favored candidate for the site of intense y-ray produc- 
tion. 

The same approach might be more difficult for MAGNs. Seven 
of the 1 1 MAGNs have a flux of F_ 9 < 40, and two of them 
(3C78/NGC1218 and PKS 0625-354) have F_ 9 £ 5. With the 
exception of NGC 1275, no significant y-ray variability has 
been d e tected so far w i th the Fermi-LAT in MAG Ns ( Abd o etall 
l2009ak iKataoka et all l2010t lAbdo et al.l l2010d) . but, at these 
low y-ray fluxes, we are unable to detect variability even when 
present. 

Significant contributions from other components can further di- 
lute and hide the flux variability. For the first time, FermZ-LAT 
has det ected y-ray emission from the lobes o f the radio galaxy 
Cen A (iFermi-LAT Collaboration et al.ll2010l) : the > 100 MeV 
emission from the lobes is a large fraction (k, 1/2) of the to- 
tal flux, with the fluxes of the northern lobe, the southern lobe 
and the core, being * 0.77 + 0.3 x 10~ 7 , 1.09 + 0.3 x 10~ 7 , and 
1.50 + 0.3 x 10~ 7 ph cirT 2 s -1 respectively. 
As discussed in Sect. 3.2.2, y-ray variability on month 
timescales has not been observed by Fermi , although flux vari- 
ability over several years is not completely ruled out. Since no 
variability has yey been clearly detected, a significant fraction of 
the y-ray flux could originate from the lobes. The y-ray spectral 
index is soft (F y = 2.52 + 0.12), implying that most of the power 
detected in the Fermi band is dominated by photons in the 100 
MeV - 10 GeV range. At these energies, the LAT PSF is unable 
to resolve the core from the lobes. As a rough estimate, we can 
rescale the y-ray flux of Cen A lobes at the proper luminosity 
distance (~ 104 Mpc) and lobe dimensions of NGC 6251. We 
use the y-ray brightest lobe of Cen A (i.e. the southern lobe) 
and assume a simple spherical geometry for the NGC 625 1 lobe 
with a diameter D - 1.9 Mpc/2 = 850 kpc, equal to half of 
the NGC 625 1 linear dimension, and thus an overestimate of the 
actual lobe dimensions. In this way, the y-ray flux rescaled to 
the NGC 6251 lobe is ,F_ 9 ss 8, comparable to the statistical 
error in the LAT flux of the source associated with NGC 625 1 
(F_ 9 = 36 + 8). There are indeed uncertainties in this estimate: 
we assume a similar SED for the lobes of the two radio galaxies 
while this cannot be assessed not even in the radio-IR band as 
the lobes of NGC 625 1 are not well sampled. For this reason, we 



discuss the effects on the nuclear SED modeling of a possible 
extended contribution in Sect. 5. 

Even though y-ray emission on the hundred parsec scale might 
also contribute, the non-thermal nuclear emission of the core 
likely dominates the s ource SED a t radio, X-ray, and y-ray wave- 
lengths. In Fig. 12 of iJones et al.l dl986l) . the core and jet-knot 
radio fluxes are compared at different frequencies: the core is 
optically thick below ~ 13.5 GHz, and becomes dominant at 
higher radio frequencies. The X-ray spectral analysis, using high 
resolution Chandra data, has shown that the bulk of the X-ray 
emission comes from the nuclear region (r < 4"), wh ile the knot 
accounts for $ 10% of the total observed X-ray flux dEvans et al.l 
2005). 



4.2. Jet orientation 

Another key question for MAGNs is the jet orientation, which 
is expected to be a large angle in objects with no extreme blazar 
characteristics. The jet inclination angle, a difficult quantity to 
observationally establish, must be measured before constraining 
the physical properties of relativistic sources. Here we consider 
a range of jet inclinations for NGC 6251 based on different ob- 
servational quantities: the jet sidedness (/), the VLBI apparent 
velocity f3„ - v a c of the jet knots, and the core dominance. If 
an intrinsic symmetry is assumed, the jet and counter-jet bright- 
ness ratio J can be expressed in terms of the jet velocity (J3 - vc) 
and orientation (0): J = [(1 +/?cos0)/(l -f3 cos 8)] (2+a \ where 
a is the radio spectral index, defined such that S (v) oc y", and 
assumed here to be a - 0.5. In a similar way the apparent trans- 



i 



verse velocity of a relativistic moving blob /?„ is related to 
and 9 via &, = ff sin /(1 - /3 cos 8) JUrrv & PadovaiH [19951 
Giovanni ni et al.l (1200 ll) found a general correlation, albeit with 
significant spread, between the core and total radio power in ra- 
dio galaxies \ogP c - 0.62 log P,„, + 7.6, where P c is the arc- 
second core radio power at 5 GHz and P, ot is the total radio 
power at 408 MHz in units of W Hz . This relation can be 
used to derive upper and l ower limits to f3 and 8 (for details, 
see Giovannin i et al.lll994l) . allowing the core density to vary 
within a factor of 2. Thus, if P c and P, ot are known, and val- 
ues or upper limits to J and j3 a are available, a narrow region of 
permitted values can be defined in the f3 versus 8 plot. This im- 
plicitly assumes that the pattern speed inferred from the apparent 
motion and the bulk motion of the emitting plasma do not differ 
(lUrrv & Pa dovani 1995, and references therein). 
In NGC 6251, t he counter jet observed in VLA maps 
(Perl ev et al.l fl984) disappears on mas scales. Hones & Wehrld 
(2002) investigated the possibility that free-free absorption by 
an ionized accretion disk could hide the receding jet, but con- 
cluded on the basis of the high electron density required that 
relativistic boosting effects better explain the observations. In 
line with their considerations, we assume that the brightness ra- 
tio of a 100:1 meas ured within 6 mas o f the core is a plausible 
lower limit to J (I Jones & Wehrld fl 9941) . The superluminal ve- 
locity of the jet has not been measured in NG C 6251, only a 
lower limit of /3 a > 1 .2 having been reported (IJones & Wehrld 
1994). The values of P c and P t0 , that we considered here to esti- 
mate the jet inclinations and velocities allowed by the core dom- 
inance correlation, were calculated fro m the 5 GHz core flux 
F c (5 GHz)=0.4 Jy JPerlevetal.il 1984 and the total 408 MHz 
flux F tof (408 MHz)=5.3 Jy {Waggett et aT] [l977h . 
The final /3-8 plot is shown in Fig. [3] The permitted values, 
shown by the gray region, are within the area delimited by 
/ > 100, j3 a > 1.2, and the core dominance relation. Only bulk 
motions larger than ~ 0.78 and inclination angles in the range of 



Migliori et al.: Jet structure in NGC6251 



00. 







1 ' i y / r r < < . 

/ / / / 












\ '•' / ' / 

\ ** f / / 


ro 




\ -**'"' / ' / 


o 




.JV' / i / 






\. v' ' / 














CO 




^^O^ ^^ 


o 




Jr**~***^^ ^t^ ^ 






*^ ^"""" , ""— — -' *^^**"^ 




^£< 






~—*^^ 






„.-■' 


p- 






o 










/? = 0sin(0) >13 






" 1-/Scos(0) 


6 




, I= [ i+^cos(e) j p> inn - 






i-/Scos(e) 


lO 




Log P c = 0.62 Log P tot + 7.6: 



20 40 

^(degree) 



60 



Fig. 3. Constraints on the angle 6 between the jet and the line of 
sight and the jet velocity B in units of the speed of light. Blue and 
black curves correspond to J - 100 and B a - 1.2, respectively. 
The red curve is from the P c versus P„„ relation. The dot-dashed 
red lines are obtained allowing the core density to vary by a fac- 
tor two. The gray area corresponds to the permitted values of B 
and 6. 



about » 10°- 40° are compatible with the observed properties 
of the NGC 6251 jet. 



5. SED modeling 

5. 1 . One-zone SSC model 

In the context of the emission from a relativistic jet, the double- 
hump shape of the nuclear SED is accounted for the nonthermal 
synchrotron and inverse Compton (IC) radiation. The radio to 
optical/UV radiation originates from the synchrotron mechanism 
by relativistic electrons, while the same electrons are supposed 
to Compton-upscatter seed-photons to X-ray and y-ray energies. 
We first model the SED of NGC 6251 with a one-zone syn- 
chrotron self-Compton (SSC) model, where the seed photons 
for Compton scattering are the local synchrotron photons them- 
selves. We assume a homo geneous spherical blob. VLBI obser- 
vations dJoneset al .11 19861) place an upper limit of S 0.20 mas 
(* 3x 10 17 cm) on the size of the nuclear region. There is no 
clear evidence of flux variability on short timescales that could 
provide more stringent limits. VLBI radio fluxes vary of close 
to two on timescales of years dEvans et al. 1 20051). correspond- 
ing to a region radius R ~ 10 18 cm ~ 1 pc. lGliozzi et al.l (12004) 
reported the detecti on of short-time (~ 10 4 s), low-amplitude X- 
ray flux variations. lEvans et al](l2005l) . analyzing XMM-Newton 
and Chandra data, concluded that variations in 2 - 10 keV fluxes 
are plausible but uncertain. 

We derive upper limits similar to the radio ones on the basis of 
our X-ray temporal analysis of XMM and Swift observations 
(Sect, 3.2.2). The assumption is that the nuclear X-ray emission 
is produced in the same region as the radio and y-ray radiation. 
As discussed in Sect. 4.2, the jet/counter-jet flux ratio, appar- 
ent velocity, and the core-dominance relation define a range of 
allowed values for the core jet speed B and relative inclination 



8. We consider as reference values for the modeling 6 = 25° 
and B - 0.91, values that fall roughly at the center of the area 
of permitted values (see Fig. |3J, The resulting Doppler factor is 
6 = 2.4. This value agrees with the constraints inferred by as- 
suming that the region is transparent to y-ray absorption via the 
yy to e + e~ proces s, 6 > 0.44 (calculated using Equation A4 in 
lAbdoetalJl2010dl Appendix). 

Looking at the broad-band SED (Fig. @}, the synchrotron 
and Compton peaks can be approximately located in the fre- 
quency intervals 10 13 — 10 14 Hz and 10 21 -10 22 Hz, respectively. 
Quillen et al. (2003) report a spectral index aos =1.1 between 
the sub-millimeter (870 /jm) and the optical band, while the IR- 
to-UV slope between 1 6,000 and 2,200 A is a ]R -nv = 1.75+0.16 
(Chiaberg e et alj|2003l) . The spectral index values in the high 
energy band, ax - 0.88 ± 0.04 in the range 2-10 keV and 
ot = 1.52 + 0.12 in the Ferm/-LAT frequency interval, are 
slightl y flatter than, but still consistent with, the sub-mm/UV a 
values. Ouille net alJ d2003l) noted that dust extinction can par- 
tially reduce the flux of the optical continuum in the core, and 
determined the spectral steepening. 

The four spectral indices give the slope below (aos, ®x) and 
above {am-uv, or) the synchrotron and Compton peaks, respec- 
tively. Therefore, as a reasonable approximation, we assume the 
same spectral indices for the synchrotron and IC curves below 
(»! = 0.88) and above («2 = 1-52) the peaks, and adopt a broken 
power-law to describe the electron energy distribution (EED) 



My) 



Ky-P' 

Ky'r" 



for y„„„ < y < y b , 
for y b < y < y max , 



(1) 



where y is the electron Lorentz factor in the range y„„„ < y < 
ymax, K is the EED normalization, p\ - 2a\ + 1 and pi - 2«2 + 1 
are the low and high energy EED spectral indices, and y b is 
the energy break of the distribution. From the ratio of the syn- 
chrotron (v s ) and IC (v f ) peak frequencies 



7b 



3v r 
4v 



1/2 



(2) 



which yields values for y b spanning the range between « 3 x 10 3 
and * 3 x 10 4 . 

The best SSC model 'fit' for the nuclear SED of NGC 6251 
is shown in Fig. |4] (left panel), and the corresponding model 
parameters (in the comoving frame) are reported in Table [2] 
Interestingly the nuclear SED can be broadly reproduced by as- 
suming a mildly relativistic motion (T = 2.4) , i.e. significantly 
lower than that required to account for several observational 
properties of blazars (i.e. the appare nt superluminal motion 



observed on milliarcsecond scales, see Kellermann et a . 2004; 



iLister et al.l [20091 and SED modeling. iTavecchio et all feoiOh . 



but with values of the source size (R - 1.2 x 10 17 cm) and 
magnetic field (B = 0.037 G) consistent with those inferred for 
BL Lac objects. Similar results have also been found when a 
single-zone synchrotron/SSC model is applied to the nuclear 
SED of other members of the MAGN sample detected by 
Fermi-LAT (see th e studies on Cen A, NGC 1275, and M87, 
Abdo et al. 2010d, 2009a b, respectively). For comparison, 
typical bulk Doppler factors and magnetic fields of BL Lac 
objec ts are 5 = 12 and B = 1 .5 x 10~ 2 G for Mrk 501 dAbdo et al.l 
2011), while for the extreme case of PKS 2155-304 during 
flaring states, detailed SSC modeling gives 5 = 100 and B = 0.1 
G(Finkeetal. 2008). 

We note that the adopted bulk motion value implies a jet 
radiation cone Q ra d = 1/T ~ 24°. Thus, the observer line of 
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sight falls roughly inside the radiation cone of the emitting 
region and NGC 6251 is aligned in the sense that it is within 
1 IT. Once again, we underline that the definition of MAGN is 
strictly related to the expectation of detecting radiation from a 
blazar-like (i.e. highly relativistic) region even in those radio 
sources whose jet physical axis is oriented at larger inclination 
angles (8 > 10°) with respect to the line of sight. 
The values of the main parameters are similar to the results o f 
the SSC model fitting of NGC 6251 in IChiaberge et al.l (120031) . 
where 8 — 18° and T - 3.2 were assumed. The particle-to- 
magnetic field energy density ratio is U' e /U' B i2 400, where 
U' e - m e c 2 n e {y) (m e is the electron mass, (y) is the average 
electron Lorentz factor) and U' B - B 2 l%n. The violation of 
the minimum energy assumption is then rather severe already 
without considering any contribution from the protons. 
In summary, the SSC model provides a good overall fit to 
the data, but requires a large departure from the equipartition 
between relativistic electrons and magnetic field an d a relatively 
slow s peed of the jet. We note incidentally that ISudou et all 
(2000) discussed a possible sub-parsec scale acceleration of the 
jet from/? ~ 0.13 ~ 0.5 pc to B ~ 0.42 at =s 1 pc. According 
to this, we should be observing the jet before it becomes 
relativistic. 

How robust are these results? The source orientation (with the 
related quantities, T and 8) is indeed a crucial assumption. A 
smaller 8, say ~ 10°, while requiring a higher bulk Lorentz 
factor, T a 15, is disfavored by radio observations: this 8 would 
imply either a large lin ear size of the sour ce, ^5.5 Mpc, or a 
significantly bent jet (Chiaber ge et al.ll2003l) . 
We note that 8 and T, as also R, affect in the same way the 
synchrotron and IC curves and will not change our results for 
the particle to magnetic field relative energy densities. Beyond 
the uncertainties related to the assumed parameters (eg 8), we 
also considered those cause by observational uncertainties and 
examined three possibilities: (i) the y-ray flux was overesti- 
mated bacause of contamination from an external contribution 
(for example the kpc jet or the lobes, as i n the case of Cen A, 
see iFermi -LAT Collab oration et al.l 1201 Ot lAbdo et al.l l2010dt 
Bottcher et al. 2008); (ii) both the y-ray and X-ray fluxes are 
overestimated, because of the limits of the X-ray and y-ray 
observatories in resolving the dimensions of the emitting region 
(R « 10 17 cm); and (iii) the flux contamination introduces 
uncertainties in the values of spectral indices . As a test, we 
considered the case where the fluxes were overestimated by at 
most a factor of 5. We assumed an intrinsically harder source, 
assigning to a\ and at the values 0.84 and 1.4 respectively, 
based on the uncertainties on the X-ray and y-ray spectral 
indices. Summarizing the results for the three cases, we have 
that B, K, and y h vary in the ranges (3.7 - 6) x 10~ 2 G, 
(1.5 - 0.25) x 10 6 cm" 3 , and (0.8 - 2) x 10 4 respectively. 
It follows that U' e IU' B may be reduced at most an order of 
magnitude but still the jet remains particle-dominated. 



5.2. Structured jet models 



Table 2. NGC 6251 nuclear SED - models. 



IChiaberge et al.l d2000l) pointed out that the luminosities of FRI 
radio sources exceed those predicted by a simple unification 
scenario, which assumes that FRI radio galaxies represent the 
de-beamed counterparts of BL Lac objects whose emission is 
produced, via SSC, by a single region. Moreover, model fitting 
of NGC 6251, as als o of other radio galaxies, do require (mo d- 
est) beaming factors (IChiaberge et al.l200lHAbdo et al.l2009allbl 



Parameters 



Model 1 
SSC 



Model 2 



Layer 



Spine 



R(cm) 


1.2 x 10 17 


8.0 xlO 16 


1.0 x 10 16 


L (cm) 


- 


2.5 x 10 16 


1.0 xlO 15 


K (cm" 3 ) 


1.5 x 10 6 


7.0 x 10 4 


7.1 x 10 4 


/min 


250 


60 


100 


Jmax 


1.4 x 10 5 


4.1 x 10 4 


2.0 x 10 4 


lb 


2.0 x 10 4 


7.0 x 10 3 


3.0 x 10 3 


P\ 


2.76 


2.76 


2.1 


P2 


4.04 


4.04 


3.0 


B(G) 


3.7 x 10~ 2 


0.7 


1.8 


r 


2.4 


2.4 


15 


6 


2.4 


2.4 


0.7 



1201 Odl) . If the condition of a homogeneous one-zone emitting 
blob is r elaxed, other possible interpretations are viable. For 
example, Ghisellin i et al.l d2005l) proposed a stratified jet struc- 
ture, with a fast narrow spine responsible for the BL Lac-like 
emission, and a slower layer, whose emission would domi- 
nate that of the fast spine at larger angles, i.e. as in FRI ra- 
dio sources. If true, there would be a strong feedback between 
the two components. The spine provides additional seed pho- 
tons via synchrotron emission that can be IC scattered by the 
relativistic electrons of the layer (and vice versa). Similar so- 
lutions, assuming a complex velocity axial/radial structure in 
jets, were also proposed by other autho rs ( Stawarz & Ostrowski 
2002; Georganopoulos & Kazanas 2003). The presence of these 
structures are also supported by the transve rsely resolved bright- 
ness p rofile observed in parsec scale jets dGiroletti et al.l 12004, 
120061) . 

In view of the new Fermi detection and additional (Spitzer and 
Planck) data in the mm-IR band, we reconsider the spine-layer 
scenar io previously proposed for NGC 625 1 in iGhisellini et al.l 
(12005b . The 'fit' to NGC 6251 nuclear SED thus obtained is 
shown in F\g^\(right panel). The geometry of the two regions is 
cylindrical (L is the cylinder length), with the spine nested inside 
the hollow layer cylinder. The spine parameters in Table [2] are 
chosen to reproduce a standard BL Lac object, with T sp j„ e = 15 
for a region of radius R spine = 10 16 cm in a magnetic field 



B 



spine 



1.8 G. The spine EED is a broken power-law with 
p\ =2.1 and p2 = 3.0 and an energy break yb, S pine — 3.0 x 10 3 . 
We keep the same spectral shape as before for the EED of the 
layer (p\ = 2.76 and p2 = 4.04), and the same 5 and 8. The 
red and black solid (dashed) curves in Fig. @] (rig ht panel) rep- 
resent the layer (spine) synchrotron and SSC contributions, re- 
spectively. The blue long-dashed curve is the layer IC emission 
of the spine synchrotron photons. The interaction between the 
spine and layer ensures the efficiency in reproducing the ob- 
served X-ray to y-ray emission and requires a low particle den- 
sity {Ki ayer = 7 x 10 4 cm 4 ), in favor of a larger magnetic field 
(Buyer = 0.7 G). On the other hand, the debeaming (6 sp j„ e ~ 0.7) 
of the spine emission at 8 = 25° hides the presence of the emis- 
sion from the fast flow. When the jet axis inclination with respect 
to the line of sight is reduced, to say 8 < 5°, the blazar-like emis- 
sion f rom the spine domi nates the SED (as also shown in Figure 
5 in Ghisellini et al. 2005, for similar parameter values). 
A critical point is that the spine-layer model assumes a flux vari- 
ability in the spine radiation similar to that observed in blazar 
sources. Rapid changes in the spine synchrotron flux should be 
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Fig. 4. Nuclear broadband SED of NGC 625 1 compiled using multi-epoch data. Radio to optical-UV data are reported in Table 
H Solid blac k triangles are the total IR flux. The solid IR points are the extrapolated non-thermal synchrotron fluxes given in 
Leipski et al. (2009). Empty optical-UV data are dereddened (see text). X-ray emission is prese nted in Sect. 3.2: black dots are for 
XMM-Netwon data a nd red from Swift satellite. Black y-ray points correspond to EGRET flux (Mukheri ee et al.ll2002l) . Red y-ray 
points are Fermi data dAbdo et al.l2010cb . Left panel: the SED is modeled with a one-zone SSC model (solid black line: synchrotron 
curve, solid red line: IC emission). Right panel: the SED is modeled with the spine-layer model illustrated in Ghisellini et al. (2005). 
Solid black and red lines reproduce the SSC emission of the layer. Dashed black and red curves are the SSC model for the spine. 
The long-dashed blue line is the IC emission of the spine synchrotron photons scattered by the layer relativistic electrons. 



reflected by those of Comptonized emission coming from the 
layer, on timescales of about t var ~ {R S pinelc)6i ayer ~ 0.02 yrs. 
In NGC 6251, no clear indication of variability in y-rays has 
been found, while the X-ray emission has varied on a shorter 
than two-year timescale. The lack of significant variability in the 
y-ray band, however, cannot by itself rule out the model. Even 
without considering the source detectability, clearly the volume 
of the layers could partly dilute the variability and of course the 
real jet structure is likely to be more complex than assumed in 
this simplified description. 

Though beyond the aim of this work, we recall that more re- 
cently two classes of models have been proposed to explain the 
high energy emission in blazars and radio galaxies. In the first 
case, beams of particles or'jets in a jet' emitting also off the 
jet axis are invoked to explain the observe d TeV flare s in both 
blazars (eg [Ghisellini et al. 2009; Giannios et"aT1l2009l) and the 
misaligned AGN M87 (Giannios et al. 2010). In the second sce- 
nario, the high-energy emission is produced by colliding shells 
with different bulk motions. In a schematic and simplified view, 
(internal) shocks are caused by the broad range of velocities 
characterizing the different "shells": the slower ones would na t- 
urally lead to emission at larger viewing angles (Dermer 2010). 



5.3. Jet power 

While from the modeling of the SED and its variability proper- 
ties it is not possible to discriminate between the SSC and spine- 
layer scenarios, these have different implications as far as the jet 
energetics is concerned as inferred here. Jet powers for electrons 
(L e ), protons (L p ), and Poynting flux (Lb) are estimated as, 



where U'. is the energy density, in the comoving frame, of elec- 
trons, U' e , cold protons, U' - m p c 2 n p (m p proton mass), and 
magnetic field, U' B . Here, we assume that there is one cold pro- 
ton per emitting electron (n p - n e ). The radiation luminosity is 



L r = nR l Y l cU' rad ~ LT 1 



(4) 



where L' is the intrinsic total (synchrotron and IC) luminosity. 
The kinetic, magnetic, and radiative powers for the two applied 
models, SSC and spine-layer, are shown in Table[3] 



Table 3. NGC 625 1 Jet Powers for the SSC and spine-layer mod- 
els. 





Model 1 


Model 2 




SSC 


Layer 


Spine 


.Mergs" 1 ) 


1.6 x 10 44 


3.9 x 10 42 


4.8 x 10 43 


L p (ergs-') 


4.9 x 10 44 


5.2 x 10 43 


2.2 x 10 44 


L B (ergs" 1 ) 


3.6 x 10 41 


2.3 x 10 43 


4.8 x 10 43 


L ki „ (ergs -1 ) 


6.5 x 10 44 


5.6 x 10 43 


2.7 x 10 44 


U (ergs" 1 ) 


2.0 x 10 43 


2.4 x 10 43 


9.0 x 10 43 



U = nR z r^cU' i 



(3) 



As noted before, the SSC model implies a strong particle 
dominance (L e * 1.6 x 10 44 erg s _1 and L p * 5 x 10 44 erg s~') 
over the magnetic field (Lb ~ 3.6 x 10 41 erg s~'). The radiative 
power is L r ^ 2 x 10 43 erg s~', about an order of magnitude 
lower than the kinetic one. In this case, the bulk of the jet 
power is conserved and goes into the formation of the large 
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radio structures. The values are in the ran ge of powers found 
by modeling the SEDs of typical BL Lacs (Celotti & Ghisellini 
2008; Ghisellini et al. 2010; Tave cchio et alj2010b . even though 
we n ote that NGC 6251 is in th e low tail of Lb values for BL 
Lacs (ICelotti & Ghisellinill2008l) . 

The SSC model - if correct - would have two interesting 
implications: (i) L e and L p are comparable to those of BL Lacs 
with similar L r , despite of the different Lorentz factor, and thus 
the jet is relatively heavier; and (ii) the low Poynting flux seems 
to exclude a magnetically confined jet. On kpc scales, the jet 
expands at variable later al velocity, and ex hibits the presence 
of re-confinement sites (Perley et al. 198 4J). Several a uthors 
dPerlev et all Il984t iMacket all 119971: (Evans et alJ 120051) have 
shown that the pressure exerted by the extended (out to 100 
kpc) halo of X-ray emitting gas around the radio source can 
account for the jet confinement on kpc scales. On the other 
hand, in the inner few arcseconds the jet is rapidly expanding 
and thermal confinement would require an X-ray luminosity 
incompatible with the observed o ne. An initial regi me of free 
expansion seems then more likely (Per lev et alJ ll984). requiring 
a re-collimation mechanism. 

In contrast, in the spine-layer model the jet is faster and 
less heavy. Thus, magnetic fields could play a fundamental 
role in its confinement. The total jet kinetic power derived 
from the spine-layer model parameters is just slightly smaller 



6. Jet power versus accretion 

The SED of NGC 6251 in the optical-UV band is dominated by 
the non-thermal emission of the jet. No sign of a big blue bump is 
identified and there is no strong evidence in the X-ray spectrum 
of a broad Fe Ka emission line, which are both considered to 
be signatures of emission by an accretion disk. The nonthermal 
flux can be used to derive an upper limit to the disk l uminosity, 
Ldisk £ 10 43 erg s" 1 dFerrarese &Fordl 19991: lHoll2009l) . which in 
units of Eddington luminosity, Lbm* is ^diskl^ -EM ~ 10 4 for an 
estimated black hole mass Mbh - 6 x 10 8 Mq dFerrarese & Ford! 
1999). The sub-Eddington regime could indicate the presence 
of a radiatively inefficient disk (here we generically refer to 
the class of models as radiatively inefficient accretion flows, 
RIAFs , see [Narayan 2002, for a review). Ghisellin i & Celottil 
(2001) demonstrated that the FRI-FRII dividing line in the ra- 
dio l uminosity vs op tical galaxy luminosity plane (proposed by 
Ledlow & Owen|[l9 96) can be expressed in terms of a critical 
accretion rate of 10~ 2 - 10~ 3 , in units of Eddington accretion, 
and suggested that this might reflect a change in the accretion 
mode from a standard o ptically thick geometrical ly thin efficient 
Shakura-Sunyaev disk ( Shakura & Sun vaevl 19731) to a RIAF. 
iFerrarese & Fordl dl999t) estimated for NGC 6251 a Bondi ac- 
cretion rate of M Bondi ~ 5 x 10~ 4 M0 yr _1 and a correspond- 
ing accretion power of P accr ,Bondi = rjM Bondi c 2 ~ 3 x 10 42 erg 



(L kin = L e + L p « 3 x 10 44 erg s 
However, in this scenario, because of the higher energy density 
in seed photons than in the SSC model, the required emitting 
particle density is low er. This leads also to quasi-equipartition 
between U' e and U' B dGhisellini et alJ 120051) . The bulk of the 
jet kinetic power is carried by the fast spine. The radiative 
dissipation (L r — 9 x 10 43 erg s _1 ) is rather high in the spine. 
The jet converts about 30% of its total power into radiation, 
and should then undergo a progressiv e deceleration, as als o 
predicted by the Compton rocket effect (Ghisellini et al. 2005), 
while we recall that NGC 6251 has an extraordinary linear 
extension of ~ 1.9 Mpc. 



5.4. Checking the jet kinetic power 

The inferences about L^,, rely on some key assumptions, the ra- 
tio of the number density of cold protons to that of electrons 
above all. A check on these estimates can be done using an al- 
ternative approach, based on the measure of the pdV work done 
by the source in forming the lobes in the hot surround ing gas 
dWillott et al.lll999UFabian et alJ l2002: Alle n et alll2006h. Here. 
we adopted the correlation reported by Wil lott et al.1 dj999) 



) than L kin of the SSC model. s ' assumin 8 an efficiency for the conversion of the accreted 



rest mass into energy of 77 = 0.1. Given the observed (non- 
thermal) luminosity, they argued that a RIAF would require an 
accretion rate 10 2 - 10 3 higher. However, as they pointed out, 
this accretion rate value is likely to be a lower limit based on 
lower limits to the pressure and density of the interstellar mat- 
ter (ISM). If we assume that the jet is completely powered by 
the accretion process, the jet kinetic power provides us with an 



> L, 



'kin 



£*«„ = 3 x 10 21 / 3/2 4£ ra-g s"\ 



(5) 



even higher lower-limi{j to the accretion power P a 

10 44 - 10 45 erg s . When we adopt typical gas densities found 
for the central regions of elliptical galax ies, rid ~ 0.1 - 0.5 crrT 3 
dDiMatteoetalJl200ll: iPellegrinil 120051) . we obtain a Bondi ac- 
cretion rate M Bon di ~ 5 X 10" 2 M© yr _1 , which gives P aC cr,Bondi ~ 

10 45 erg s -1 . If this is the case, a radiatively inefficient disk can in 
principle provide the required jet power, though the mechanism 
that converts the accretion into the jet power is very efficient. A 
study of nearby X-ray luminous elliptical galaxies have found 
that P aC cr,Bondi and the jet power, inferred from the work done to 
expand the cavities observed in the surrounding X-ray gas, have 
similar values. This suggests that a significant fraction of the 
mass entering the Bond i accretion r adius eventually goes into 
the relativistic jet (Allen et al. 2006). Interestingly, it has been 
propo sed for blazars (C elotti & Gh isellini 2008; Ghisellin i et al.1 
2010, and references therein) that the jet power maintains a lin- 
ear dependence on M/M£dd with a transition to a quadratic de- 
pendence at low accretion rates, as also expected in RIAF solu- 
tions dNaravan & Yilll995UMahadevanlll997h . 



where Z451 [erg s _1 Hz sr _1 ] is the monochromatic radio power 
at 151 MHz. In the revised formula considered here, a factor / 
takes into account possible systematic underestimates intrinsic 
to the technique. lHardcastle et al.1 (120071) estimated / to be in the 
range between 10 and 20 for a sample of FRI and FRII sources. 
Given the 151 MHz luminosity of NGC 6251, L151 = 2.46xl0 31 
erg s _1 Hz -1 sr , the kinetic power, setting / respectively equal 
to 10 and 20, is in the range 7.5 x 10 43 -2.2 x 10 44 erg s _1 , just 
slightly below the value of L^,„ found for the two models, thus 
in agreement with the above estimates. 



7. Summary 

We have presented a study of the broad-band nuclear SED of 
the radio galaxy NGC 6251, one of the MAGNs detected by 
Femw-LAT. In agreement with previous studies, the nuclear 
SED is dominated by non-thermal emission related to the sub-pc 
jet. A nuclear origin appears to be the most likely hypothesis for 



5 Indeed, this lower limit relies on the assumption that there is a pro- 
ton component, which seems however supported by independent jet 
power estimates (see Sect. 5.4). In addition, we assume that the pro- 
ton component is not due to the jet entrainment of intergalactic matter. 
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the y-ray emission. This is also supported by the results of our 
X-ray analysis of archival XMM-Newton, Chandra, and Swift 
observations. Both models, SSC and spine-layer, adopted to 
reproduce the observed SED, provide a good overall fit. In the 
following, we summarize the main results of the SED modeling. 

A single-zone SSC model assumes a slow and heavy jet. 
The particle power dominates over the magnetic field power 
by about three orders of magnitude. This seems to rule out 
a magnetically accelerated and confined jet. The X-ray halo 
surrounding the jet can account for its confinement on kpc scales 
but the core region remains poorly constrained. A low bulk 
Lorentz factor, when a one-zone SSC model is adop ted, is also 
shared w ith other MAGNs detected by Fermi-LK ^ (lAbdo et al.l 
l2009allbT: JKataoka et alj|20!ol lAbdo et alj|2010dl) . This suggests 
the presence of at least a second site of y-ray emission in 
addition to the one we observe w hen the jet is pointing directly 
toward us (Chiaber ge et alj|20 00). A significant reduction of the 
jet inclination angle, which would allow higher bulk Lorentz 
factors, is unsupported by the radio observations. These results 
should not be significantly affected by an overestimate of the 
y-ray flux of a factor of a few (« 5). 

A structured jet with a fast inner component and a slower layer 
is a possible explanati on. In the spine-layer model proposed in 
iGhisellini et al.1 (120051) . the jet is relatively light and could be 
magnetically confined. The fast inner component carries the 
bulk of the jet power. On the other hand, the strong radiative 
dissipation (L r ~ 0.3 L^ m ) is at odds with the Mpc-length of the 
jet and the flux variability predicted by the model in the high 
energy band has not yet been observed. 

The derived jet powers, as 10 44 - 10 45 erg s -1 , are model- 
dependent and rely on some important assumptions. However, 
similar values wer e derived using the 151 MHz luminosity 
(IWillottetal.lll999h . 

The disk component appears to be completely hidden by 
the jet emission. The sub-Eddington luminosity regime 
(Ldisk/LEdd ^ 1CT 4 ) could be related to the presence of a RIAF. 
Support to this hypothesis is provided by the estimates of the 
Bondi accretion rate, albeit with significant uncertainties in 
the ISM physical parameters. As has been found for blazar 
sources, the mechanism channeling the accretion power into the 
jet seems to work in a more efficient way {Lunl^Usk i£ 10 and 

'-'kin ~ U.i — lr accr)- 
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Table 4. NGC 6251 nuclear SED - data. 
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Frequency 
(Hz) 



Flux 

(Jy) 



Flux error 

(Jy) 



Reference 



Angular 
Resolution 



3.26 x 
4.09 x 
6.09 x 

1.40 x 

1.41 x 
1.67 x 
2.30 x 
2.30 x 
2.70 x 
4.75 x 
5.00 x 
5.00 x 
8.40 x 

1.05 x 
1.07 x 
1.50 x 
1.50 x 

1.54 x 
2.28 X 
3.00 x 
3.30 x 
4.00 x 

4.07 x 
4.30 x 

6.08 x 
7.00 x 
9.35 x 
1.00 x 
1.43 x 
2.17 x 
3.45 x 
1.00 x 
1.30 x 
2.00 x 
2.00 x 
1.00 x 

1.42 x 
1.85 x 

2.55 x 
3.69 x 
5.40 x 
5.45 x 
8.73 x 
1.24 x 
1.52 x 
3.00 x 
5.37 x 
9.58 x 

5.06 x 
1.66 x 



10 s 

10 s 

10 s 

10 9 

10 9 

10 l > 

10 9 

10 9 

10 l > 

10 9 

10 9 

10 l > 

10 9 

10'° 

10 10 

10'° 

10'° 

10 10 

10'° 

10'° 

10 10 

10'° 

10'° 

10'° 

10'° 

10'° 

10'° 

10" 

10" 

10" 

10" 

10' 3 

10' 3 

10' 3 

10' 3 

10 13 

10 14 

10' 4 

10' 4 

10' 4 

10' 4 

10' 4 

10' 4 

10' 5 

10 22 

10 22 

10 22 

10 22 

10 23 

10 24 



0.70 
<0.30 

0.60 

0.50 

0.36 

0.28 

0.25 

0.65 

0.72 

0.60 

0.36 

0.13 

0.72 

0.55 

0.90 

0.34 

0.66 

0.85 

1.30 

1.34 

1.50 

1.35 

1.50 

0.50 

1.20 

0.84 

0.70 

0.82 

0.70 

0.47 

0.42 
4.97 x 10~ 2 
9.00 x 10- 2 
2.79 x 10- 2 
0.81 x 10~ 2 
1.64 xl0~ 2 
1.08 x 10~ 3 
5.45 x 10"* 
3.14 xlO"* 
2.00 x 10"* 
1.34x10"* 
1.21 x 10~ 4 

3.30 x 10- 5 
1.77 x 10" 5 

< 1.62 x 10~ 10 
4.10 x lO" 11 
1.85 x 10-" 
1.83x10-" 
1.23 x 10~ 12 

1.31 xlO~ 12 



0.10 

0.10 
0.10 



0.02 
0.03 



0.01 
0.05 



0.01 
0.04 
0.04 
0.06 
0.07 
0.13 
0.08 
0.02 
0.10 
0.15 
0.30 
0.09 
0.08 
0.07 
0.01 

2.00 x 10~ 



0.11 x 10~ 3 
0.54 x 10- 4 
0.31 x 10- 4 
0.20 x 10~ 4 
0.13 x 10~ 4 
0.12 x 10~ 4 
0.33 x 10- 5 
0.18 x 10~ 5 

3.07 x 10-" 
1.12 x 10-" 
0.84 x 10 -11 



7 
6 
7 
6 

6 
7 
6 

7 
7 
7 
8 
9 
10 
9 
9 
9 
11 
11 
11 
11 
11 
11 
11 
11 
12 
12 
12 
12 
12 
12 



WENSS (n) (55") 

3.7' 

WENSS (o) (56") 

NVSS (fl) 

VLBI (<3 mas) 

VLBI (0.6 mas) 

3.7" 

VLA {fl) 

VSOP (0.5 mas) 
VLA (a) 



VLBA (0.5 mas) 

VLA (fl) 

0.65" 

Wmap (0.88°) 

Planck (32.65') 

0.66° 

Planck (27') 

0.51° 

VLA (a) 

0.35° 

Planck (13.01') 

0.22° 

Planck (9.94') 

Planck (7.04') 

Planck (4.66') 

23" 

4-11" Total emiss 

0.5-2' 

4-11" Total emiss 

4-11" Sync. emiss. (w 

4-11" Sync, emiss. *' 

0.2" (c) 

0.2" (c) 

0.2" w 

0.2" w 

0.2" (c) 
02 "(O 

0.2" <c) 

0.2" (c) 

0.97°x0.74° 

0.97°x0.74° 

0.97°x0.74° 

0.97°x0.74° 

0.97°x0.74° 

0.97°x0.74° 



Notes. (a) Archival data. {b) Extrapolated IR synchrotron flux (Leipski et al. 2009, see text) . w Optical-UV non de-reddened flux. 



References. (1) Wagget tetafl ( fi977h: (2) Hones et alj (fl986l): (3) [C ohen & Readhead ( 1979); (4) Sudou & Taniguchi (2000); (5) Mac k et al 
( 1997); (6) Wri ght et al.U2009h : ( 7)[Ma nck Collabor ation et alJd201lh ; (8)lQuillen et aill l2003l) ; (9)[Leipski et al. (2009); (10) K napp et al.ldl99i 
(ll) IChiaberge et al.1 (12003); (12) Mukherie e et al.ld2002h . 
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